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ABSTRACT 


A  method  of  determining  two-dimensional  pressure  distributions  on 
arbitrary  foils  is  explained.  The  development  is  based  on  an  approximate 
potential  theory  suggested  by  Moriya  (1941)  which  *s  empirically  modified  in  a 
manner  suggested  by  Pinkerton  (1936)  to  give  an  arbitrary  lift  for  a  set 
incidence  and  at  the  same  time  satisfy  the  Kutta  condition.  Interpola¬ 
tion  functions  for  the  ordinates  are  used  to  reduce  the  calculation:,  to 
a  straightforward  numerical  procedure  which  is  easily  programmed  for 
machine  calculations.  (FORTRAN  statements  for  such  a  computer 
program  are  included.)  The  results  are  those  of  Riegels  (1943). 

Comparisons  are  made  with  othe»  theoretical  methods  for  poten¬ 
tial  flow  and  with  experimental  results.  Good  agreement  with  both 
calculated  and  measured  pressure  distributions  is  found  when  the  lift 
coefficients  ore  matched.  The  assumption  that  cavitation  occurs 
when  the  local  pressure  falls  to  the  vapor  pressure  is  not  upheld  for 
the  cases  considered. 

ADMlNISTPv^TIVF  INFORMATION 

This  report  is  essentially  a  duplication  of  o  Master's  Thesis  submitted  to  the 
Graduate  School  of  Cornell  University  in  May  1965.  The  work  was  performed  at  the 
Taylor  Model  Basin  under  Bureau  of  Ships  Subproject  S-R009  01  01 . 
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INTRODUCTION 


To  obtain  on  accurate  estimate  of  the  actual  pressure  distribution  on  two- 
dimensionci  fo-'n>s,  we  have  two  approaches:  (1)  to  solve  the  cumbersome  Navier- 
Stokes  equations  governing  viscous  flow  or  (2)  to  use  a  potential  theory  with  empirical 
modifications  to  approximate  real  fluid  effects.  The  first  approach  may  be  simplified 
to  solving  the  boundary- layer  equations  for  a  two-dimensional  curved  surface,  but 
the  simplification  is  nominal  and  the  task  remains  formidable  in  application.^'^ 
Moreover,  as  yet,  this  method  is  not  sufficiently  advanced  to  give  accurate  results 
unless  experimental  boundary-layer  data  are  known.®  The  second  possibility  is  based  on 
potential  theory  which  means  simpler  develcpment  and  shorter  computation  time.  As 
normally  used,  this  method  makes  use  of  the  experimental  lift  which  is  either  available 
or  can  be  estimated  for  conventional  foils. 

The  first  satisfactory  method  of  using  the  experimental  lift  in  potential  theory 
was  developed  by  Pinkerton  in  1936.'^  He  made  use  of  the  Betz  (1915)  observation^ 
that  employing  the  experimentally  detcmined  circulation  in  place  of  the  theoretical 
value  gives  a  pressure  distribution  thot  agrees  well  with  measured  results  over  most  of 
the  chord.  However,  merely  inserting  the  experimentally  determined  circulation 
into  potential  theory  will  not,  in  general,  satisfy  the  Kutta  condition  that  the  flow 
leave  the  trailing  edge  smoothly.  Pinkerton"^  was  able  to  retain  the  measured  lift 
and  also  satisfy  the  smooth-flow  condition  by  introducing  an  empirical  modification 
in  the  profile  shape.  He  applied  the  modification  to  a  rigorous  potential  theory 

^  References  ere  listed  on  page  72. 
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formulated  by  Theodorsen  in  1932.'  Although  Theodorsen's  theory  is  easier  to  apply 
than  the  bcundoiy  layer  equations,  the  calculations  are  still  lengthy, 

Riegels"^  removed  this  last  complication  in  1943.  He  adopted  Pinkerton's  idea 
of  an  empirical  distortion  in  the  flow  field  to  an  approximate  potential  theory 
developed  by  Moriya  in  1933.'  Riegels  made  a  further  simplification  by  linearizing 
the  equations  on  a  small  parameter  Introduced  in  the  modification.  In  addition,  he 
used  interpolation  functions  for  the  ordinates  which  reduce  the  calculations  to  a 
straightforward  numerical  procedure,  requiring  only  offsets  at  fixed  fractions  of  the 
chord.  The  linearized  numerical  procedure  of  Riegels  has  been  in  use  for  some  time* 
at  the  David  Taylor  Model  Basin  where  it  was  programmed^  for  the  computers  at  the 
Applied  Mathematics  Laboratory.  In  addition  to  giving  the  chordwise  pressure 
distribution,  the  results  of  this  program  can  be  used  to  predict  cavitation  inception 
if  the  assumption  is  made  that  the  cavitation  starts  when  the  local  pressure  falls 
to  the  vapor  pressure  of  the  surrounding  liquid.  Unfortunately,  output  from  this 
program  consists  of  the  pressure  distribution  at  only  the  input  points  so  there  is  no 
assurance  that  the  minimurr.  is  obtained.  During  the  analysis  necessary  to  correct 
this  omission,  !’■  became  obvious  that  the  Riegels  derivation  ’.'as  quite  obscure.  A 
detailed  derivation  of  the  theory  was  accordingly  undertaken,  and  the  results  of  the 
investigation  are  outlined  in  this  paper.  The  derivation  closely  follows  later  work 
of  Moriya^"  (1941)  which  is  an  approximate  conformol  transformation  of  the  circle 

*ln  1955,  a  Model  Basin  Memorandum  (Aero  28)  described  a  method  of  calculating 
pressure  distributions  over  profiles  of  arbitrary  shapes.  This  memorandum  included  a 
translation  of  the  Riegels  paper  (Reference  7). 
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to  on  airfoil  profile  and  gives  the  same  equations  for  the  velocity  distribution  as  his 
earlier  work. 

The  present  contribution  to  the  development  is  o  unified  presentation  of  the 
theory,  simplification  and  extension  of  the  numerical  method,  and  extensive  com¬ 
parisons  with  experiment  and  tieory.  In  the  comparisons  with  theory,  some  new 
results  of  considerable  simplicity  are  derived  from  an  exact  conformal  transformation 
of  the  circular  cylinder  (see  ADpendix  A).  Airfoil  geometry  is  discussed  in  Appendix  B. 

The  FORTRAN  statements  of  a  new  computer  program  for  calculating  the  pressure 
distribution  on  an  arbitrary  profile  are  given  in  the  paper  together  with  a  description 
of  the  program  and  input  instructions.  The  program  statements  are  listed  in  Appendix  C. 
This  program  will  accept  foil  ordinates  at  arbitrary  stations  and  interpolate  between 
them  to  obtain  the  ordinates  at  the  required  stations. 

POTENTIAL  PRESSURE  DISTRIBUTION  FROM 
APPROXIMATE  CONFORMAL  MAPPING 

SUMMARY  OF  DERIVATION 

Using  the  method  of  conformal  transformation,  we  can  eosily  transform  the  known 
flow  about  a  simple  body  to  that  about  an  arbitrary  profile.  As  is  customary,  we  take 
the  circular  cylinder  for  the  simple  body  since  both  its  geometry  and  complex  potential 
are  simple  in  form. 

This  development  retains  the  two  simplifying  approximations  made  by  Moriya^'’: 

(1)  a  group  of  transformation  coefficients  is  ignored  in  the  calculations  and  (2)  the 
circulation  is  determined  by  placing  a  stagnation  point  at  the  trailing  edge  of  this 
approximate  transformation  rather  than  the  exact  one.  For  symmetrical  foils,  the  second 
approximation  is  exact. 
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COEFFICIENTS  OF  THE  TRANSFORAAATION 


A  difficulty  v/ith  all  transformations  of  the  circle  into  an  airfoil  is  that  of  matching 
the  coefficients  of  the  transformation  with  those  of  the  oirfoil.  Except  for  speciol 
cases,  iteration  is  required.  In  the  v/ork  that  follov/s,  v/e  neglect  the  iteration  and 
instead  evaluate  the  coefficients  v/ith  a  first  approximation  to  the  transformation. 

The  most  general  mapping  function  which  maintains  the  type  (but  not  necessarily 
magnitude  and  direction)  of  the  flow  or  large  distances  from  the  body  is  given  by 


y>  cs 

z  =  c  A  4-C  -  zc 

-la  o  I  n 


t 


=  X  -r  ly 


n 


where  z  is  the  complex  coordinate  in  the  profile  plane. 


C  =  A  +  i  B  , 
n  n  n  ' 

n  =  -1,0,1,2, . . 


C  is  the  complex  variable  in  the  circle  plane,  and 

^  =  re'^  ,  as  shown  in  the  figure  below  which  also  shows  the  orientation  of  the 


flow  velocity. 
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Substituting  C  =  ae‘'',  v.e  obtain  >'  one  y  on  the  profile  05 


1 

< 

II 

X 

CQSC-  B_'i 

sifxj;  - 

"  (  A^  cos 

3^  sin  nv') 

'2' 

Y  =  B 

-1 

cost?-  A 

-1 

sifv:)  ~ 

3  - 
0 

T  (B  cos  np  - 
1  n 

A  sin  rvp) 
n 

'3‘ 

where  the  capital  letter  Y  Is  used  to  denote  points  on  the  profile.  To  ensure  rhot 
these  expressions  actually  give  the  ordinotes  of  the  profile,  the  coefficients  of  the 
series  must  match  a  similar  expansion  for  the  profile  coordinates.  Usually  v/e  have 
0  ^  x<  1  and  Y  =  Y{x).  In  Appendix  A,  severe!  foils  ore  derived  from  exact  theory. 
For  foils  v.'hich  are  slender  In  the  horizontal  direction,  it  is  shown  that  >;  is 

X  =  1  2  (1  -  cosci  -  0(t/0 

where  T  ond  fore  the  thickness  ratio  and  camber  ratio,  respectively.  For  such 
slender  foils,  a  slight  change  in  the  x  volue  makes  an  even  smaller  change  in  the 
ordinate  (except  in  regions  of  large  curvature,  e.g.,  the  nose).  With  these  facts  in 
mind,  we  moke  the  approximation*  that 

X  =  1/2(1  +  cosrp)  T4~‘' 

(note  that  this  approximation  reduces  ^  to  a  parameter)  and  Y  =  Y(0)  where  Y(o) 
is  understood  to  be  an  implicit  Fourier  expansion,  i.e., 

CO 

Y((d)  =  a  +  S  (a  cosri©+b  sin  n(p)  fSol 

o  ,  n  n 


*Alternatively,  this  could  be  the  starting  point  for  on  iteration  leading  to  the  exact 
transformation.  However,  the  results  show  this  approximation  to  be  sufficiently  accurate 
for  the  shapes  considered  here  (see  Appendix  A). 
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Mafching  coefficients  with  Equation  ^31,  we  obtain 


Bo  =  =  1/2-  Yd(p 

o 

2-  2- 

B  ^B  =a  =1/-  r  Ycoi(pd‘p,A  -A  =b  =  1/-  J  Y  sir!0d(£> 

-111  o  -111  "o 

2- 

B  =  a  =  1/-  Y  cos  n-rd^p  ,  “^  =  =  1/-  ^  Y  sin  rvpdrp 


n  >  2 

Also,  from  Equation  4~,  we  knov/  that  A  =  i  '?  ’ience  me  coefficienrs  are 

determinen  witr  trs  exceoricr-  /-  or.u  L  . .  ^NOf^ng  that  the  infintre  se''ies  in 

Eauctionr  "2'  one  ”2  ori=^  '’  •oner  conjugares  ano  tnof  the  A  .  and  &  ,  terms  ore  not, 

v/e  can  immccj  J'S  '  write  an  exoression  ro'  derermining  these  HTi,  riv»'^*cvc  y 

tne  resulrinc  exoressions  would  involve  A  one  E-  .  Tne  deoendence  on  rnese 
'  c  c 

coefficients  may  be  eliminated  by  considering  differentials  of  Equations  ^2l  and  ^31. 
Therefore,  expressions  which  involve  only  A  and  B  ,  may  be  obtained  by  constructing 

•*»  I  —»  • 

the  quantities:  ^ 

T  (rp)  =  dx/dcp  -  ]/2-  f  Y'(t)  cot  -  t)/'2  dt 

o 


2^ 

S  (<p)  =  dY/d(p-r  1/2-^  f  x'(t)  cot('p-t)/2  dt  r5i 

o 

0 

where  f  indicates  thor  the  Cauchy  Principle  Value  of  the  improper  integral  is  to  be 
taken  ond  the  primes  denote  differentiation  with  respect  to  the  argument.  Moriya's 
notation'/  differs  from  ours  since  he  used  the  Fourier  development  which  is  more 
cumbersome  tnon  the  integral  and  different  al  form  piesented  heie. 


SubsHtution  of  the  tronsformafion  equations  (Equations  r21  and  [3"l)  into 
Equation  fSl  yields 

T  (73)  =  -2  TA  ^  sin  (p-»  B  ^  cos  cpl 

S  (rp)  =  2  [A  ^  costp-  B  ^  sin(p^  ^5a^ 

With  the  given  airfoil  ordinates,  the  expressions  for  T  and  S  are 

2TT 

T  {(p)  “  -1/2  sin  ^  -  1/2^  f  Y'  (t)  cot  (<p  -  t)/2  dt 

o 

S  (^p)  =  l/2cos<o+ (©)  r6] 

Equating  Equations  fSal  and  [6]  for  CD=  0,  we  obtain 

=  1/4  0  +2Y'(0)1 
2rT 

B_i  =-1/4tt  J  Y'(t)  cot  t/2  dt  m 

o 

Ncic:  Y  (Oj  =  (dY/d(p)  _  =-l/2sin(p(dY/dx)  is  zero,  for  example, 

-  <p-U  TE 

with  profiles  having  Siiarp  trailing  edges.  Since  not  all  profiles  are  of  this  form, 
Y'(0)  will  be  retained  in  the  development. 

Although  the  Fourier  conjugates  method  is  only  one  possible  way  of  obtaining 
the  coefficients,  it  is  shown  later  that  the  group  of  coefficients  T(<p)  and  S('p)  do 
appear  in  the  expression  for  the  velocity  distribution.  Even  more  important,  we  have 
found  that  other  methods  of  obtaining  only  A  ^  end  B  ^  do  not  give  results  which 
compare  as  favorably  with  exact  results  as  do  the  T(<p)  and  S(fp)  expressions. 
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LIFT  COEFFICIENT 


Arguments  based  on  physical  concepts  lead  to  the  conclusion  that  there  is  no 
flow  around  the  trailing  edge  (Kutta  condition).  Thus,  the  numerical  value  of  the 
circulation  is  just  enough  to  place  a  stagnation  point  at  the  location  on  the  circle 
which  transforms  to  the  trailing  edge  of  the  profile.  In  the  exact  transformation, 
the  value  of  7)  which  makes  x  a  maximum  describes  this  locotion.  In  the  approxima¬ 
tion,  the  maximum  value  of  x  is  at  <p=  0. 

The  velocity  in  the  circle  plane  is  found  by  differentiating  the  complex  potential 
F(C)  with  respect  to  For  a  circular  cylinder  inserted  in  a  uniform  stream  such 
that  the  stream  makes  an  angle  with  the  real  axis,  the  complex  potential  is 

F(0  =  V  +  \  (T /2-)  In  ItM 

where  P  is  the  circulation. 

With  0  as  the  independent  variable,  the  velocity  on  the  circular  cylinder 
C=ae'  becomes 

dF/dCI,  +i(r/2Ta)  e"'®  [81 

C=ae  -* 

Applying  the  Kutta  condition,  we  put  dF/dC  =  0  at  C  =  a,  which  sets 

r  =  A-nQ  V  sin  (a^)  [91 

Before  the  circulation  can  be  converted  to  a  lift  coefficient,  the  effect 
of  the  mapping  function  on  the  free  stream  must  be  investigated.  In  the  circle 
plane,  the  free  stream  velocity  is  given  as 

dF/dC  I - -Ve"'°t 

^  —no 


9 


and  in  the  profile  plane  as 


dF/dz  I  = 

Z 

We  define 

dF/dz  J- 

and  let 

r  -  'J 
C  —  \e 

-1 

But 

^-1  = 

Thus 

K  =  'Va 

-lay 


-la 


Z  —  CO 


-1 


-1 


X  =  orctan 


^-1 


Equating  the  two  descriptions  for  the  freestream,  we  find 

V  =  KU/a 


=  a-x 


[101 


ml 


With  the  above  relations,  the  lift  coefficient  is  found  to  be 

pur  ! - 

"  l/2PU^c  “  2tt  V  (4  +  (4  sin  (a  -  x) 


[121 


where  c  is  the  chord  and  p  the  fluid  mass  density.  Hence,  x  is  the  angle  of  zero 
lift  in  potential  flow. 

VELOCITY  AND  PRESSURE  DISTRIBUTION 


The  velocity  on  the  profile  q  is  found  by  evaluating  the  expression 


!  dC  \a  cto/  '  ^ 


10 


Performing  the  operations,  we  find 


Substituting  Equation  ^11]  and  expanding  the  right-hand  side,  we  obtain 


nai 


+(A  cosp- B  sin^  -  A  )  sin  a  I  ri41 
-1  -1  -1 

The  combination  of  transformation  coefficients  enclosed  by  parentheses  is  the 
previously  defined  S  (tp)  and  T  (<d),  Equation  fSl.  Hence,  the  expression  for  the 
nondimensional  velocity  can  be  written 

_q_  _  I  [  T  (<P)  ~  T  (0)  1  cos  [  $  (<o)  -  S  (0)  2  sin  g  | 

U  -JT/4  sin^  (0  +  Y'{^)  2 


s  I  C  {<p)  cos  a  +  D  (0)  Sin  a  ! 
where  T  {(D)  and  S  (rp)  are  given  by  Equations  r6’lo 

Since  this  method  is  for  steady  irrotational  flow,  Bernoulli's  equation 
holds  and  gives  the  pressure  distribution  in  coefficient  form  as 


C 

P 


P  -  P 

oo 

1/2PU2 


1  - 


ns] 
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EMPIRICAL  MODIFICATION  TO  ACCOUNT  FOR  VISCOUS  EFFECTS 


In  Reference  4,  Pinkerton  shows  that  on  empirical  modification  to  Theodorsen's 
exoct  potential  theory^  gives  results  which  agree  reasonably  v/e!l  with  experimental 
measurements.  In  this  chapter,  we  will  show  hov/  Pinkertc:’'?  idea  can  be  extended 
to  the  approximate  potential  theory  just  outlined. 

Theodorsen’  noted  thal  conventional  profiles  are  similar  in  shape  and  that 
applicotion  of  the  Jourkowski  tronsformotion  to  an  arbitrary  profile  will  yield  a 
nearly  circular  shape  of  variable  radius  1/4  e^  and  polar  angle  cd.  The  nearly 
circular  shape  can  then  be  trar'sformed  into  a  circle  of  constant  radius 
1/4  e  and  polar  angle  0  +  e  {to).  With  our  notation  and  the  circulation  as  yet 
unspecified,  Theodorsen's  expression  for  the  velocity  on  the  profile  can  be  put  in 
the  form: 


where 


=  f  (^)  I  sin  (a-rp-  e)  - 


rt-e 


toy 


H(p)  = 


1  + 


and 

X  =  1/2  l^cosh  +  cosh  cosipj 
Y  =  1/2  sinh  0  sin  7? 

(Note  that  this  tp  is  different  from  that  used  in  the  previous  section.)  The  expressions 
for  X  and  Y  are  sufficient  to  determine  il>  and  and  from  these  we  can  compute 
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IS 


2n 

e  {(d)  =  i/2-  p  y(r)  cot 
o 


(f)  -i  e{o)  -  [t  +  e(0 


1  - 


dc 
dt  J 


dt 


2^ 

=  1/2-  ;  C(t) 


1  + 


dc 

dt 


dt 


Pinkerton's  modification'^  is  to  add  an  arbitrary  function  to  e  (o).  Arguing  that  the 
modification  is  a  result  of  viscous  effects  ond  ihor  the  cumulative  effect  of  the 


viscous  forces  are  toward  the  trailing  edge^  he  takes  the  addition  to  c  as 


a  = 


R 


where  is  to  be  determined  from  the  requirement  that  the  lift  be  the  experimental 
value  ano  R  is  some  constant  >  1/2  (if  R  <  1/2  the  velocity  goes  to  infinity  at  the 
nose).  Pinkerton  considered  only  the  cose  R  =  1 .  Since  C  is  the  Fourier  conjugote 
of  e,  a  change  in  c  would  also  cause  a  change  in  b.  However,  Pinkerton  has 
founo  that  the  effect  of  this  change  in  is  negligible  in  f  (o). 

For  slender  foils,  both  c  ond  are  small  and  we  can  approximate 

X  ss  1/2  (1  +  cos 
Y  ^  1/2  J/)sinO 


\  d^  dy)  j  ^ 

Thus  the  new  expression  for  velocity  is 


H+' 


_dc 

d0 


1  + 


da 


d(p  J 


■3-  =  ffe)  (1  +  — 

U  \  d(p 


sin 


(n[  _  g ')  _  a  _  ^  1  _  sin  f  a-  X-  ‘a) 


where  X  =  c  (C)  when  the  Kutta  condition  is  satisfied.  Comparing  this  expression 
with  Equation  [13],  we  can  immediately  write  a  similar  modification  to  the 
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approximate  theory  given  previously: 


( 

q  ^  2K'^1  +  do  /I  sin  (a  -  a"  -  x)  -  ssn  (; 

TT  "  _ ~~  ^ 


-  X-  la'  I 


or 


3-  =  1  + 


±l' 

do  , 


T  (o)  cos  (a-  cl)  +  S  {'p\  sin  (a  -  aO 
-  [T  (0)  cos  (a  -  t3.)  +  S  (0)  sin  (a  -  ^.a)  ^ 


dY^ 


sm  «p  ^ 

4  &p 


[16] 


where 


da^  La  .  _  cos  V 


(0\  R-1 


Except  for  notjtion,  this  is  the  solution  obtained  by  Riegels.^ 

Implied  in  t!ie  "jbove  equation  is  that  ♦he  lift  coefficient  is 

(a-  X  -  Ha)  [171 

We  seo  at  once  that  ta  is  a  fictitious  decrease  in  ihe  geometrical  angle  of  attack 
to  give  the  desired  lift.  This  equation  con  be  rearranged  to  give  6a  in  terms 
of  an  experimental  lift  coeffic'ent 

(exp) 

La  =  a-x-arcsin  _  [^g] 

in  ^(4  A_^)=  +  (4B_^)‘* 

With  known  values  of  the  lift  coefficient,  this  expression  determines  La.  Since 
we  rarely  have  specific  test  data,  a  representative  actual  lift  curve  slope 
coefficient  and  angle  of  zero  lift  can  be  used  to  find  the  lift.  Assuming  these 
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values  ore  known,  we  can  determine  the  lift  coefficient  from  the  expression 


l) 

^  'exp  s 


19' 


where  V  is  equal  to  dC.  /d3.  j  /  2-  and  o:  is  the  experimental  angle  of 

L  SXp  O0 

zero  lift.  (tjoiXi  «  are  discussed  further  later  in  the  paper.) 

In  order  to  more  clearly  show  the  effect  of  the  velocity  dist'-ibution  will 
be  linearized  on  to. .  This  will  be  put  ‘n  the  form  of  the  potential  velocity  plus 
a  term  multiplied  by  lo  .  The  I inearizotion  is  performed  by  replacing  sin  x' 
by  a  ^  and  cos  by  1 .  Ignoring  all  smoll  quantities  multiplied  by  x' ,  we  find* 


where  E  {<■/))  =  (1/2)  .  . .  -  x  D  ('/') 

Vl/4  sin^O+  Y  '  {'pT 


and 


Ixr) 

oa  =  0 


is  the  potential  velocity.  In  addition  to  showing  the  effect 


of  to.  on  the  potential  velocity,  this  expression  is  easier  for  hand  calculations 
with  the  potential  theory  known. 


NUMERICAL  METHOD 

TRIGONOMETRIC  SERIES  TO  REPRESENT  TrlE  ORDINATES 


To  use  Equation  ’^16"',  we  need  an  analytical  expression  for  the  ordinates.  For 
tabulated  offsets,  this  means  either  constructing  a  Fourier  series  or  curve  fitting  a 
polynomial  to  the  ordinates.  As  will  be  shown  in  this  section,  a  trigonometric 


*With  R  =  1,  the  linearized  expression  for  velocity  would  be  the  some  as  that 
obtained  by  Riegels^  in  1943  if  he  had  recognized  that  C  (o)  sin  a  in  our  E  ('p) 
term  is  small  in  the  normal  incidence  angle. 
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series  may  be  fitted  analytically  to  ordinates  given  at  certain  specific  stations  along 
the  chord.  The  slope  and  cotangent  integral  then  become  a  summation  of  the  product 
of  constant  coefficients  and  the  ordinates  at  these  fixed  stations. 

The  formulas  used  to  approximate  the  ordinates  can  be  developed  by  either  the 
least -square-error  approach’^  or  by  requiring  the  series  to  pass  through  the  given 
points.’'^  Both  methods  are  equivalent  so  we  will  outline  the  more  direct  method  of 
curve  fitting  through  the  given  points. 

Values  for  the  ordinates  are  required  at  2N  1  equally  spaced  values  of  o 
between  Q  and  2" in  Equation  .4  ' ,  i.e.,  at  the  fixed  stations 


_  1  +  cos  ^ 
^m  ~  - 


1  +  co: 


,  0  <  m  <  2N 


^21 


Althcugh  we  have  2N  +  1  values  of  the  ordinates,  only  2N  are  independent  since 

the  periodicity  requires  Y(0)  =  Y(2-"),  or  the  mean  value  if  there  is  a  discontinuity 

at  the  trailing  edge.  Noting  that  sin  is  zero  at  all  m,  we  see  that  Equation  [3a^ 

should  be  truncated  to  the  2N  terms 

N-1 

Y('p)  =  aQ+  r  ( 

1 

Setting  this  expression  equal  to  the  2N  independent  ordinates,  we  obtain  a  system  of 
linear  equations  to  be  solved  for  the  coefficients: 

N-1 

These  equations  are  solved  by  multiplying  both  sides  of  the  above  equation  by  the 
multiple  of  the  desired  coefficient  and  summing  from  zero  to  2N  -  1 .  Using  this 


a  cos  nrp 


+ 


b„sin 


x\tp)  +  a  cos  Nco 
'  N 
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approach  we  obtain  the  constants 


I 


1 


2N-1 


a  =  — —  _  Y 

o  2N  0 


1  2N-1 

°n°Tr  5 


1  2N-1 

a. ,  = -  I  Y  cos  m- 

N  2N  0  "" 


b  =  i- 

n  N 


2N-1 


0 


mn- 

'TT 


(Note  that  for  n  these  expressions  are  identical  v/ith  results  of  evaluating  Integral 
expressions  for  the  Fourier  :onstants  by  the  trapezoidal  rule  with  equal  angular  spacing 
as  above.)  With  Y(0)  =  0  =  Y(-),  the  expression  for  the  ordinates.  Equation  „3a  , 
becomes 


'22' 


where 


N-1  ^  N-1  N-1 

1  m=l  n=l 


23‘ 


N  1  N-1 

Y  (co)  =  r  a„  cos  n^J  =  “  "  (Y  -  Y^. .  ) 
even  ■  n  ^  M  '  m  2N-m 

0  m=1 


N-1 

n=l 


cos  n  'r  cos  n  v 
’  m  ' 


1  cos  N^pj^  cos  Np 


'24' 


ond  Y^  and  Y2|Ng_^  ore  the  upper  and  lower  surface  ordinates,  respectively,  meosured 

m~ 

from  the  nose-tail  line  at  fp=  =  -jq-  . 

The  odd  part  of  the  series  arises  from  thickness  and  the  even  part  from  camber. 

It  should  be  noted  that  the  odd  function  given  above  assumes  a  rounded  leading  and 
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d(Y  +  Y2n_  J 

trailing  edge  ond  the  even  function  assumes  finite  slope  - - -  at  the 

dx 

leading  and  trailing  edges. 

EVALUATION  AT  SPECIFIED  POINTS 

The  troublesome  expressions  in  T  {'p)  and  S  (rp)  (Equation  [6])  are  the  singular 
corangent Integral  and  the  slope,  respectively.  Using  the  numerical  expressions  for 
Y  (rp)/  we  can  easily  perform  the  necessary  operations.  The  expressions  v/ill  involve 
a  double  summation  as  in  the  Y  ifp)  expression.  Since  the  summation  involving 
the  trigonometric  terms  Is  to  only  a  finite  limit,  the  sum  may  be  calculated 
analytically.*  At  a  point  corresponding  to  an  input  point  {(p=<p^  ^ 


m),  the  cotangent  integral  reduces  to 

}"  N-1 

1/2tt  Y'(t)  cot “2  =  V2N  Z 

o  m=1 

Y 

2N-,n 


1  -  (~1) 


m+ 1; 


Y 


m 


m-u 

L_cos-|q- 


r-1 


cosH^  -1 


(for  m  =  Uf  the  multiple  of  Y^  is  easily  shown  to  equal  N/2). 


Similarly,  the  slope  (in  the  fp  plane)  reduces  to 


dY  \  ,  ^  \  J  +  m+  V  r. ,  (m  -  j/)  TT 

<i«  l'P=  'Pj,  m=l  L  2N 


(m  +  u)  TT 

2N-m''°^  2n' 


(for  m  =  i/,  the  multiple  of  Y^  is  zero). 

The  numerical  expressions  of  Riegels^  may  be  added  together  to  obtain  these  results. 
However,  these  equations  require  computation  of  a  N  +  1  by  N  ~  1  array  of  coefficients 

*See,  for  example.  References  13  and  14. 
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this  is  a  lengthy  calculation  either  by  hand  or  on  a  high-speed  computer.  Further 
simplification  is  possible  if  we  define  >  =  m  -  i/  in  the  first  term  of  summotion  and 
X  =  m  +  y  in  the  second.  The  equations  then  reduce  to 


O  -t 

l/2’-|Y'(t)cot-2^ 


^  V 


1 

2N 


N-] 


X=1 


Y  ~  Y 


1  -  (-1) 


cos 


>- 

N 


-] 


dY\ 


=  1/2 


A=1 


Y  -  Y 

Vx 


cot 


X- 

2N 


'25" 


:26'’ 


where  ws  understand  that  Y  =  Y^. ,  and  Y^k  =  Y  . 

-m  2N-m  2N+m  m 

These  expressions  are  now  quite  simple  and  can  easily  be  set  up  for  hand  or 
computer  calculations  with  a  minimum  of  calculations.  It  is  necessary  that  only  the 
ordinates  be  cyclically  rearranged  before  each  computation. 

The  expressions  for  T  {'p)  and  S  {<p)  ai  (p  =  tp  are  now 


T^-T(<pJ  =  -l/2 


N-1 


N 


) 


1  -(-]) 


smS  +NY,-i-  -  (Y  +Y  . 

N  1-cosi^ 


X=1 


"N  J 


S  hS(c.^  =  1/2 


cos  -jq-  + 


N-1 


X=1 


(  ^-x-Vx  W 


and  the  modified  velocity  at'p=(0  is 

V 


27' 


^28] 


3. 

U 


T  cos  (a  ~  a  0  +  3  sin  ( a  -  n  0 

v  y 


“'"To  cos  (a-  -’a)  +S^  sin  (a  -  ‘.a)] 


X 


y 


\/ 1/4  sin^^  —  + 

'  N  d(D  / 


'P 


-o.y- 

1  -  -y  sm-jq- 


y 


1 6a^ 
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EVALUATION  AT  AN  ARBITRARY  POINT 


Ai  menUoned  in  the  Introduction,  the  motivation  for  this  research  is  to  obtain 
the  minimum  piessure.  Unfortunately,  there  is  no  guarantee  that  the  pressure 
distribution  computed  at  the  fixed  points  will  yield  the  minimum  value.  If  we  note 
that  only  when  the  minimum  is  near  the  nose  does  the  pressure  distribution  hove  o 
steep  gradient,  then  our  problem  can  be  reduced  to  obtaining  more  data  lust  near 
the  nose. 

The  procedure  outlined  for  obtaining  the  slope  and  cotangent  integral  at 
specified  points  can  be  followed  for  an  arbitrary  tp  value.  The  expressions  obtained 
are  more  complicated  but  still  give  a  single  summation  of  the  ordinates  at  the  input 
points  and  trigonometric  expressions  w-.ich  can  be  evaluated  once  and  for  all  at  any 
desired  point.  Although  the  results  appear  to  be  new,  the  derivation  is  not  sufficiently 
involved  to  warrant  presentation  and  wc  give  only  the  final  results: 


2-  j 

1/2- P  Y'(t)  cot-^-dt 


N-1 

=  E  (Y  -Y 

-1  ^ 
rci —  I 


2N-m 


) 


rn 


(- 1 )  cos  Np  -  I 


:N 


sin^^'^m  sin^ 


(cosrp-  COS(p^)^ 


/  .  M  5in<p 

(-1)  Sin  Nrp  ^  m 


COSfp  -  COS  'p 


m  -I 


N-1 

m=l 


m 


(-1)  COS  N^p  -  1 

2N 


1  -  cos^^cosc/? 

(coS'^  -  cosO 
'  m' 


m+  1 

(-1)  sin  N<p 


sin  <D 


cos  (D  -  cos  <p 


[291 
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Figure  1  —  Chordwise  Stations  =  1/2  (1  +  cos  m-/lS) 


TABLE  1 


\  aIuos  of  ^  1/2(1  cos  mrr  18) 


m 

0.  deg 

* 

0 

0 

1.0 

1  > 

0.S92404 

2 

20 

0.969846 

3 

30 

0.933013 

< 

40 

0.883022 

5 

SO 

0.821394 

6 

80 

0.75 

7 

70 

0.671010 

8 

80 

0.586824 

9 

90 

0.5 

10 

100 

0  413176 

n 

110 

0.328990 

ii 

120 

0  25 

130 

0  178606 

M 

UO 

0  116978 

IS 

ISC 

0  066987 

le 

160 

0  030154 

17 

170 

0  007596 

'  IS 

180 

0 

*36-in  *m 


\ 

» 


I 


jv  N-1 

—  =  ^  Y  -  Y  ) 

*>  ra=l  "■ 


s!n  M/)  .  5'"'(’mSin.p 


2N 


(cos0  -  COS7D^)=^ 


I  M  s:n<P 

+  (~*)  Htp  .  m 


cos^  -  cos<P 


m 


N-1 

r 

m=  1 


+  v  (Y  +Y..,  ) 
m  2N-m 


/  ...  ^  -  cos<P  cos^j 

(-1)  sm  N^)  .  m  ^ 


2N 


(cos  //)  -  cos  f/5 
'  m  ' 


m  +  1 

(-1)  cos  N<y)  . 


sm  tp 


cos  (p  -  cos  tp 


m 
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In  order  to  obtain  a  better  estimate  of  the  minimum  pressure,  the  velocity  is 
evaluoted  at  additional  points  near  the  nose  using  these  expressions.  Choice  of 
these  additional  points  depends  upon  the  choice  of  N.  in  the  pressure  program,  N 
will  be  18,  so  thato  increases  in  increments  of  10  degrees;  see  Figure  1  and  Table  1 . 
The  intermediate  points  selected  are  listed  in  Table  2.  They  are  each  integer 
degrees  between  the  nose  {<o=^)  and  the  input  point  on  each  side  of  it 
{(0  =  tt  ±  )  and  three  additional  points  2  1/2  degrees  apart  between  this  point 


(0=Tr  ±  -H—  )  and  its  neighbor  (rp  =  tt  ±~—  ). 
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COMPARISONS  WITH  ANALYTICAL  RESULTS 


The  numerical  expressions  for  the  ordinates  are  obtained  by  forcing  the 
expression  to  pass  through  specified  points.  Thus  there  is  no  assurance  that  the 
computations  v/ill  be  reasonable  between  the  points.  Also,  we  have  assumed  that  a 
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finer  spacing  of  computed  data  near  the  nose  is  sufficient  to  determine  the  minimum 
pressure.  We  cannot  check  these  problems  in  a  general  sense,  but  we  can  make 
comparisons  in  specific  cases. 

In  order  to  check  the  cccuracy  of  the  numerical  solution,  the  computations 
were  performed  both  analytically  and  numerically  for  the  linearized  symmetrical 
Joukowski  profile^^ 

X  =  1/2(1  +  cos 


added  to  the  parabolic  camberline^' 


(sin  o  -  1/2  sin  2d) 


Y  =  fsin^O 
c 

where  T  is  the  thickness  ratio  and  f  is  the  camber  ratio.  The  analytical  and  numerical 

solutions  were  accurate  to  at  least  four  significant  figures  when  the  input  was  specified 

to  five  significant  figures.  This  accuracy  Is  more  than  adequate  for  any  practical  use. 

In  order  to  make  a  comparison  with  the  minimuiTi  pressure  computed  analytically 

and  numerically  at  discrete  points,  the  results  for  a  symmetrica!  ellipse  were  used. 

For  this  simple  shape,  the  location  and  the  value  of  the  minimum  pressure  are 

easily  obtained  analytically.  From  the  exact  potential  velocity  distribution 

(see  Appendix  A),  we  obtain  the  location  of  maximum  velocity  =  0  | 

\dfp  / 


from  the  following  equation 


tan  a  ~ 


sin  75(1  -  f  )-  tan'p 
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where  r  is  the  thickness  ratio.  Porficular  o  values  can  be  selected  and  the  angle  of 
incidence  for  maximum  velocity  calculated  .from  this  equotion.  A  graphical 
comparison  of  both  the  location  and  value  of  the  maximum  velocity  as  computed 
analytically  and  numerically  at  the  fixed  points  is  given  in  Figure  2  for  thickness 
ratios  of  0.05  and  0.1.  This  comparison  shows  that  even  though  the  location  in 
the  nose  region  may  be  off  somewhat,  the  velocity  is  calculated  with  sufficient  accuracy. 

COMPUTER  PROGRAM 

A  computer  program  has  been  written  for  calculating  the  pressure  distribution 
by  the  numerical  method  outlined  in  the  previous  section.  The  FORTRAN  statements 
of  the  source  program  are  listed  in  Appendix  C.  These  statements  are  for  Applied 

Mathematics  Laboratory  Problem  840-041 -F. 

The  order  of  the  computer  operations  is  as  follows:  First,  the  constant  coefficients 

in  Equations  [27]  through  [30]  and  other  constants  are  calculated  and  stored.  Then  the 

input  data  are  read  and  the  summed  products  of  the  ordinates  and  the  stored  constants 

obtained.  Finally,  for  each  angle  of  attack,  the  pressure  distribution  at  the  input 

points  and  at  the  intermediate  points  is  determined  from  Equation  [16^,  (actually  Equation 

[16a]  of  the  numerical  method)  with  a,'  =  6ax{i.e,,  R  =  l).  Since  the  lost  operation 

performed  is  for  tne  angle  of  attack,  variations  in  angle  of  attack  take  least  machine  time. 

The  data  needed  to  compute  the  pressure  distribution  consist  of  the  angles,  in 

degrees,  for  which  the  pressure  distribution  is  desired;  estimates  of  the  angle  of  zero 

dC  » 

lift,*  in  degrees,  and  an  average  lift-curve  slope  ccefficient*  v  =  — /2tt 

da  / exp 

*Reference  15  finds  that  77  is  primarily  a  function  of  the  thickness  form  and  that 
is  primarily  a  function  of  the  meanline.  Tabulated  values  for  many  sections  can  ® 
be  found  in  Reference  16.  As  shown  by  the  data  in  Reference  16,  both  77  and  aQ 
are  functions  of  the  Reynolds  number.  ® 
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LOCATION  OF  MAXIMUM  VELOCITY 
»  -  0.  DEGHEES 


TABLE  2 


for  use  in  Equation  ^19'';  and,  of  course,  the  cirfoil  ordinates,  in  fractions  of  the 
chord.  The  nondirr.ensional  ordinates  should  be  occurate  to  about  one-tenth  of  1 
percent  of  the  maximum  ordinate  to  give  a  srrooth  pressure  curve.  (Accuracy  to  four 
decimal  places  is  usually  sufficient.) 

The  stations  for  whicl  the  offsets  ore  -  squired  in  the  numerical  method  are 
given  by  Equation  *^21  ],  shown  in  Figure  1, and  tabulated  in  Table  1.  Since  these 
stations  are  often  inconvenient  to  use  with  tabulated  data,  provision  fo'  input  at 
arbitrary  x  values  is  also  provided.  The  input  at  the  required  stations  is  obtained 
from  a  third-degree  polynomial  fitted  through  four  ordinates,  two  on  each  side  of 
the  required  station.  The  interpolation  is  performed  between  angular  stations  since 
the  ordinates  do  not  then  have  on  infinite  slope  ot  the  leading  edge  (see,  e.g.. 
Figure  10).  The  arbitrary  stations  should  be  as  near  the  required  values  as  possible. 
This  is  especially  important  at  the  leading  and  trailing  edges. 

As  mentioned,  several  input  options  are  provided.  As  with  all  machine 
computations,  input  must  be  given  in  a  nonvarying  monner;  hence  the  input 
descriptions  must  be  rigidly  followed. 

Option  1;  Ordinates  at  Required  Stations 

The  shortest  method  of  data  input  is  to  specify  ordinates  at  the  required  stations. 
This  type  of  input  naturally  results  in  the  shortest  running  time  since  no  interpolation 
for  the  ordinates  is  required.  The  necessary  order  of  input  is  shown  in  Figure  3a 
where  each  horizontal  block  represents  a  different  card  starting  from  Column  1 , 
However,  not  all  the  caros  indicated  for  the  angles  of  attack  and  the  ordinates 
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should  be  filled;  only  os  many  boxes  os  are  needed  should  be  used,  (See  Jhe  sompie 
input.  Figure  3b.)  When  several  profiles  ore  given  for  the  same  angles  of  attack, 
it  is  important  that  only  similar  types  of  profiles  be  grouped  together  (i.e.,  either 
symmetric  thickness  forms  or  nonsymmetric  cambered  foils,  not  both).  VVhen  changing 
types  of  foils,  o  nev/  set  of  input,  starting  with  the  initial  control  card  for  the  ongles, 
is  necessary.  Additional  members  of  a  group  of  profiles  require  only  the  cards 
indicated  under  PROFILE  INFORMATION  in  Figure  3a. 

Option  2:  Ordinates  at  Arbitrary  Stations 

As  mentioned,  input  may  also  be  given  at  arbitrary  stations.  Two  sub¬ 
options  are  considered  in  this  case. 

(a)  Points  on  Profile  » 

The  first  of  these  sub-options  is  to  specify  points  (x,Y)  around  the  foil 
in  the  direction  of  increasing  m  (see  Figure  1).  The  upper  surface  trailing  edge 
ordinate  must  be  the  first  point  specified  and  the  lov/er  surface  trailing  edge  ordinate 
must  be  the  last  point  given.*  The  nose  (see  Equation  B2^  in  Appendix  must 
also  be  given  before  the  input  points.  In  addition,  the  nose  point  should  be  given 
again,  in  order,  midway  in  the  listing  of  coordinates.  The  program  shrinks  and 
rotates  the  coordinatss  to  put  the  nose  at  (0,0)  since  this  orientation  has  been  found 
to  be  the  most  occurate.  (This  point  is  further  discussed  in  Appendix  B.)  The  rotation 
angle  is  added  to  each  of  the  input  angles  so  that  a  is  measured  from  the  original 
reference  line.  A  maximum  of  53  points  is  permitted. 

*Fven  tliough  there  ore  provisions  for  inserting  ordinates  ct  the  trailing  edge,  aid  though 
they  are  printed  out  in  the  pressure  program,  no  use  is  made  of  them  when  computing  the 
pressure  distribution.  They  are  used  in  the  interpolation  and  hence  shv^jld  not  be  omitted 
when  using  the  arbitrary  input  option. 
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Figure  3a  —  Format  for  Data  Input  to  Computer  Program  at  Reo.uued  Stations 
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Figure  3b  -  Sample  Input  at  Required  Stations 

Figure  3  —  Data  Input  to  Computer  at  Required  Stations 
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(b)  Values  of  Thickness,  Camber,  and  Ccmberltne  Slope 

The  second  sub-option  is  to  give  values  of  the  station,  the  thickness,  the 

camber,  and  the  camberline  slope  from  trailing  edge  to  leading  edge.  These  quantities 

are  combined  in  the  NACA  manner’  (thickness  applied  perpendicular  to  the  comber) 

to  obtain  the  coordinates  of  the  foil  (see  page  1 13  of  Reference  15  and  also  Equation 

_B1~  in  Appendix  B  of  this  report).  The  coordinates  are  rotated  and  shrunk  as  above. 

If  the  surface  is  to  be  computed  by  Y  -  Y  :::  Y  ,  the  camberline  slope  should  be  entered 

c  t 

as  zero  in  the  input  and  the  nose  radius  is  not  needed.  Input  at  27  stations  Is  permitted. 

The  format  for  data  input  in  both  cases  above  is  shov/n  in  Figure  4a,  and 
sample  input  In  Figure  4b. 

Other  Options 

Although  not  shown  in  Figures  3a  ond  4a,  if  a  1  Is  placed  in  Column  16  of  the 
control  cord  just  before  the  PROFILE  INFORMATION  cards,  values  of  the  lift 
coefficient  (as  many  as  the  angles  and  the  same  format)  may  be  given  in  place  of  the 
angle  of  zero  lift  and  lift-curve  slope  coefficient.  (See  the  sample  input.  Figure  4b.) 
Note  that  lift  coefficients  must  be  given  for  each  profile  grouped  under  this  control 
card.  Also,  if  a  1  is  placed  in  Column  l2ofthe  initial  control  card,  only  values 
of  the  angle  of  attack,  lift  coefficient,  minimum  pressure  coefficient,  maximum 
velocity,  X  location,  and  integrated  moment  coefficient  about  the  line  x  =  0.25 
are  printed  out  instead  of  the  complete  pressure  distribution  around  the  foil.  There 
is  a  considerable  saving  in  paper  and  computer  time  when  using  this  option. 
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Figure  4  —  Data  Input  to  Computer  at  Arbitrary  Stations 
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Figure  4a  —  Format  for  Data  Input  to  Computer  at  Arbitrary  Stations 
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Figure  4b  —  Sample  Input  at  Arbitrary  Stations 
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Sample  output  for  the  RAE  101  (10-percent  thickness)  is  shown  in  Figures  5  and  6, 
This  is  the  data  plotted  in  Figure  15.  The  first  two  pages  of  output  (Figure  5)  for  each 


set  of  ordinates  are  titled  Profile  Constants.  The  first  page  includes  the  stations  (column 

headed  X)  and  the  ordinates  (column  heoded  Y).  The  column  headed  C  is  the  velocity 

distribution  on  the  profile  at  zero  angle  of  attack  and  the  column  headed  D  is  the 

velocity  distribution  due  to  the  angle  of  attack  (Equation  '^15'').  The  column  headed 

E  is  a  factor  necessary  to  find  the  ruodified^  linearized  pressure  distribution  (see 

Equation  '^20]),  These  three  columns  permit  rapid  evaluation  by  hand  of  the 

velocity  at  other  than  the  input  angles  without  rerunning  the  data.  (At  the  bottom 

of  the  page  is  an  expression  which  shows  how  to  combine  Columns  C,  D,  and  E 

dY  dY 

to  obtain  the  linearized  pressure  distribution.)  The  last  column  is  -  -  1/2  sinf^  r — 

jr-  cr®  dx 

The  theoretical  slope  factor  - - /  2’^  I  is  printed  out  at  the  bottom  of  the 

da  Cj_  =  0 

page  as  well  as  the  theoretical  angle  of  zero  lift. 

Additional  information  may  be  obtained  from  the  data  on  the  Profile  Constants 

page  of  output.  The  angle  of  attack  at  which  a  stagnation  point  lies  ot  the  nose,* 

(q)  ft  “  0/  calculated  for  potential  flow  from  Equation  ri5 

TT  /  Z.  a  -  0 

C(n)  cosa+  t)(’^)  sin  a  =  0 


i.e. , 


a  =  arctan  - 


C(tt) 

D(rr) 


For  0,  the  nose  stagnation  angle  may  be  similarly  obtained  from  the  linearized 


*This  angle  is  not  the  same  as  Theodorsen's  ideal  angle  of  attack/  The  ideal  angle 
(for  thick  foils)  is  defined  as  the  angle  of  which  a  stagnation  point  lies  at  the  forward 
end  of  the  camberline.  As  noted  by  Theodorsen,  this  angle  is  of  limited  practical 
importance. 
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guro  T)  —  T\pu-al  Output  from  ('omputor,  Profili'  Cori'-tant 
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Figure  6  —  Typical  Output  from  Computer,  Pressure  Distribution 
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velocity  given  by  Equation  '^20'' : 

C(’- )  cos  0.  )  sin  c:  *  laE(-)  =0 

Usually  two  iterations  are  necessory  to  solve  this  equation. 
Also,  if  the  nose  radius  is  desired,  it  moy  be  found  from 


dffi 


1  p 


(-) 


if  the  nose  radius  is  known,  this  relation  provides  a  check  of  the  computed  slope. 

Two  pages  titled  "Pressure  Distribution"  are  printed  (Figure  6)  for  each  Input 

angle.  The  first  is  the  pressure  distribution  at  the  specified  points  and  the  second 

is  the  pressure  distribution  at  the  intermediate  points  near  the  nose.  Constants  printed 

out  at  the  top  of  the  page  are  for  the  input  conditions.  "Delta"  stands  for  ia  , 

calculated  from  Equation  18*',  and  is  given  in  radians.  A  measure  of  the  accuracy 

dY 

of  the  input  data  is  the  smoothness  of  the  curve  since  small  errors  In  the  ordinates 


are  magnified  when  taking  the  slope. 


At  the  bottom  of  the  page  with  the  pressure  distribution  at  the  input  points 


(Figure  6)  are  integrated  values  for: 

2^ 


C 

c 


=  -  1/2^  C  sin  (pd'p 
o  P 


2-^ 

=  -#  CpY'(o)d(p 


2- 

C..  I  =^1/2  0  C  sin'p  (x-.25)  d'P 

x=.o  O  “ 

2- 

C  I  =-|  C  YY'MdO 
y  =  0  °  P 


The  coefficient  of  force  normal 
to  the  chord! ine,  (positive  upward) 

The  coefficient  of  force  parallel 
to  the  chord! ine  (positive  in  the 
positive  x-direction) 

Moment  about  the  line  x  =  0.25 
(clockwise  positive) 
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Moment  about  the  line  Y  =  0 
(clockwise  positive) 


The  lift  coefficient  is  given  by  C,  =  C. ,  cos  a  -  C  sina^  and  the  moment  coefficients 

^  '  L  N  c  ' 

about  the  point  (0.25,  0)  are  given  by  C  (0.25,0)  =  C  I  +  C  I 

x=0.25  ^  Y  =  0 

The  integrals  are  evaluated  v/ith  the  trapezoidal  rule  ct  equal  angular  intervals. 

The  numerically  Integrated  values  should  be  reasonably  accurate,  except  since 
there  is  not  a  smoothing  factor  (e.g.,  sin  cd,  Y)  under  the  integral  sign.  In  fact, 
the  pressure  program  gives  a  nonzero  chord  force  for  foils  at  zero  incidence  . 

(Potential  theory  gives  zero  drag.) 

In  all  coses  considered  so  far,  the  integroted  and  set  lift  coefficients  have  agreed 
within  2  percent, 

COMPARISON  OF  POTENTIAL  THEORY  WITH  OTHER  METHODS 
Here,  comparisons  will  be  made  with  other  methods  for  calculating  two-dimensional 
pressure  distributions  in  potential  flow.  Substituting  Equation  [6]  into  Equation  riSl, 
we  find  for  the  potential  velocity  on  the  profile 


q. 

U 


_ 1 _ 


COS  a 


1/2  sin^ 


2rr 


-l/2Tr^  Y'  (t)  ^cot  ^  -  cot  t/^t 


1 


+  sin  a 


1  -  cos  tp 

- 2 -  -  Y'  {(P)  +  Y'  (0) 


This  is  equivalent  to  the  solution  obtained  by  Moriya.®'^® 


[311 
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Converting  this  expression  to  x  dependence,  we  obtoin 


q. 

U 


1  Y/(0-Y/(0 

-  I  2-  P  — - _  d^ 

o  X  -  4 


where  0  <  x  •"  1,  and  Y  and  Y  ore  the  upper  and  lov/er  surface  ordinates,  respective 
-  -  U  L 

measured  from  the  nose-tail  line,  and  primes  denote  differentiation  with  respect  to 
the  argument. 

Riegels  and  Wittich’"^  obtained  this  same  equation  by  noting  that  linear  airfoil 
theory,  including  the  thickness  correction  term  in  sin  a,  giv'cs  a  velocity  on  the 
chordline.  They  obtained  the  velocity  on  the  profile  from  the  linearized  velocity 
q(x,0)  by  equating  integrands  of  two  expressions  for  the  circulotion: 


i.e 


/ 


r  =  .9q(x,Y)  ds  =  q(x/0)  dx 

,  wx  q(^/0)  =  q(x,o) 

q(x, Y)  s;  - “  — ■  - 

ds/dx  1+  Y  (x)^ 


since  the  singularities  enclosed  by  both  paths  are  the  same.  This  neat  trick  can 
also  be  used  to  improve  three-dimensional  linear  theory.^®' 

uquurion  ’’32]  cleorly  shows  the  cr^' :s-coupl ing  effects  of  thickness  and  camber 
in  the  velocity  distribution.  If  the  thickness  and  camber  are  odded  and  subtracted  to 


obtain  the  ordinates,  then  the  only  coupling  at  zero  incidence  is  the  term 


■yl  +  (x)"  which  is  negligible  except  near  the  ncse.  Hence  the  NACA 

method’  of  combining  velocity  increments  computed  separately  for  thickness  and 
camber  can  be  expected  to  be  accurate  everywhere  except  near  the  nose.  Comparative 
calculations  have  siiown  that  using  the  NACA  tabulated  velocity  increments^"  is 
sufficiently  accurate  for  most  sections  but  that  there  is  simply  not  enough  data  near 
the  nose  even  for  moderate  incidence. 

In  References  13  and  14,  Weber  proposes  a  method  of  calculating  two-dimensional 
pressure  distributions  v/hich  is  based  on  other  work  of  Riegels.  Her  expressions  for 


C  (<p)  (Equation  is  the  same  as  ours,  bi^t 


D  (^)  H  0  (x) 


1 

1  f* 

1 +—P 
-  o 


or  D  ('p)  she  obtains 

Yp'n 

""irfny 


d^ 

x-C 


V  l  +  Y  '  (x)  ' 


where  Y^  =  1/2  (Y^  -  Y^). 

This  expression  gives  results  which  are  closer  to  exact  potential  theory  than  does 
the  method  presented  earlier  in  this  paper.  However  the  derivation  is  somewhat 
loose  when  camber  is  included  end  the  method  does  not  easily  allow  the  insertion 
of  the  experimental  lift  while  satisfying  the  Kutto  condition. 

The  data  tabulated  in  Reference  15  will  be  used  to  compare  the  approximate 
potentioi  solution,  Equation  ^32/  with  the  exact  solution  for  the  more  usual  foils  in 
use  today.  Theoretical  velocity  contributions  for  a  large  number  of  arbitrary  thickness 
forms,  calculated  numerically  from  Theodorsen's  exact  method  for  potential  flow,'' 
are  tabulated  there.  Pressure  distributions  at  zero  angle  of  attack  calculated  by  that 
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method  and  by  the  method  of  this  report  are  corripcred  for  three  foils — the  NACA 
16-009,*  the  NACA  0010,  and  the  NACA  65A.-006.  As  shown  in  Figures  7  ond  8, 
the  two  methods  are  in  close  agreement  for  both  the  NACA  16-009  and  the  NACA 
0010.  Hov/ever,  for  the  65  form  (Figure  9)  the  method  of  this  report  predicts  a  small 
pressure  peak  near  the  nose.  Plotting  the  ordinates  at  the  angular  variable 
n  =  arc  cos  (2  x  -  1)  instead  of  the  usual  x,  we  find  a  hump  near  the  nose  (Figure  10). 
This  hump  causes  the  pressure  peak  since  fairing  it  out  was  found  to  give  a  smooth 
pressure  curve.  The  numerical  method  of  Hess  and  Smith^  predicts  a  similar 
peok  on  this  foil  (see  Figure  42  of  Reference  18).  They  note  the  important  point 
that  experiments  fail  to  shov/  the  peak.  The  experimental  discrepancy  is  thought 
by  the  authors  to  be  the  result  of  inaccurate  machining  near  the  nose. 

To  explain  the  theoretical  differences,  note  the  discussion  on  page  3  of  Referencv»  19 
in  connection  with  the  design  of  the  NACA  16  series: 

“The  Theodorsen  method  as  ordinarily  used  for  calculating  the 
pressure  distributions  about  airfoils,  was  not  sufficiently  accurate  near 
the  leading  edge  for  prediction  of  the  local  pressure  gradients." 

Because  of  this  deficiency,  Theodorsen's  numerical  evaluation  of  a  nonlinear  integral 
equation  was  not  used  for  the  design  of  the  NACA  6  series  foils,  but  was  replaced  by 
another  numerical  procedure  employing  equations  similar  to  the  interpolation  functions 

*The  equations 

Yy  =  T  r. 98883-/^-  .23702  x  -  .04398  x""  -  .55762  x®  1  ,  0  <  x  <  .5 

Yj  =  T  F. 01  +  2.325  (1-x)  -  3.42  (l-x)2+  1.46  (l-x)3  1  ,.5  <x<^  1 

give  a  good  fit  to  the  tabulated  ordinates  for  the  NACA  16  form.  The  nondimensional 
nose  radius  is  0.4889  r^,  or  0.003960  for  the  original  9-percent  thick  foil. 
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for  the  ordinates  (Equations  r23l  and  [24^)  used  in  this  report.^/  For  these  foils, 
it  is  believed  that  the  complicated  behavior  of  the  "design  parameters"  near  the  nose 
(see  Figure  4  on  page  9  of  Reference  19)  is  not  adequately  represented  by  such  a 
finite  term  approximation  for  the  Fourier  series  unless  a  large  number  of  terms  is  used. 
Note,  in  cc.itrast,  the  smoothness  of  the  curves  Y  (o);  see,  e.g..  Figure  10.  Accurate 
integration  of  an  irregular  curve  would  depend  upon  the  number  of  points  taken. 
Theodorsen®  recommended  five  equal  divisions  between  0  and  with  his  method. 
Reference  15  recommends  40  with  the  improved  method,  surely  more  than  enough. 

No  indication  of  the  number  used  in  the  design  process  is  given  .‘n  Reference  19, 
but  it  is  believed  to  have  been  insufficient. 

So  far,  comparisons  have  been  made  only  for  thickness  distributions  at  zero 
angle  of  attack.  For  cambered  foils  at  an  angle  of  attack,  a  potentiol  theory 
(to.  -  0)  comparison  can  be  made  or  Aamoy  be  adjusted  to  give  the  same  lift.  Both 
comparisons  are  shown  in  Figure  1 1  for  the  441 2  at  a  =  6.4  degrees.  (Pinkerton‘S 
gives  the  potential  lift,  computed  numerically  from  exact  theory,  as  =  6.915 
sin  (a  -  a^),  aQ  =  -0.0706  radians.)  It  is  easily  seen  that  the  computed  pressure 
distribution  with  the  exact  lift  (C|^  =  1 .254)  is  not  in  agreement  with  the  results 
of  the  exact  potential  theory.  It  should  be  noted  that  the  results  of  the  exact  theory 
were  taken  from  the  small  figure  on  page  62  of  Reference  15.  Even  allowing  for  errors 
in  the  transfer  of  data,  the  agreement  is  only  fair  at  the  same  lifts.  Moreover,  the 
compjted  minimum  pressure  is  about  10  percent  lovver  than  the  "exact, "  which 
follows  the  trend  in  Appendix  A.  Integration  of  the  curves  in  Figure  11  gives 
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Figure  10  —  Ordinate  \ersu?-  Angular  N  ar^able.  Showing  Dump  n«‘ar  'he  Leading  Kdge,  NA.C'A  65A008 


•> 


=  1 .21 A  for  the  ''exact"  pressure  distribution  and  =  1 ,242  for  the  opproximote 
case  (icc  '  0),  compared  to  the  set  of  1.254.  The  difference  between  the  set 
and  integrated  value  for  the  pressure  program  (approximately  1  percent)  is  thought  to 
be  caused  by  the  approximations  in  the  theory  and  possibly  by  inaccuracies  in  the 
numerical  method.  The  disagreement  between  the  "exact"  results  and  the  integrated 
value  (approximately  4  percent!  is  obviously  CJ*<sec?  ir  part  by  the  small  size  of  the 
original  curve  of  exact  ''or.-  'vhtch  the  data  v/ere  taken.  Some  of  the 

discrepancy  could  al'o  be  a  result  of  Pinkerton's  use  of  approximate  expressions 
for  e  and  That  these-  are  approximations  is  apparent  if  one  compares  the  expressions 
on  page  13  of  Reference  4  with  those  in  Reference  6  and  in  the  section  on 
empirical  modification  given  in  the  present  paper. 

COMPARISON  WITH  EXPERIMENTAL  RESULTS 


CO/'^f^ARISON  WITH  MEASURED  PRESSURE  DISTRIBUTIONS 


There  are  few  measured  pressure  distributions  tabulated  at  a  sufficient 
number  of  points  to  accurately  define  the  pressure  curve  near  the  nose.  The  most 
comprehensive  tests  are  those  of  Pinkerton^”^  on  the  NACA  4412  section.  Unhappily, 
Pinkerton's  measurements  are  actually  for  three-dimensional  flow  over  a  rectangular 
foil  with  a  30-inch  span  and  a  5-inch  chord.  The  equivalent  two-dimensional  flow 
is  found  by  subtracting  the  theoretically  calculated  induced  angle  of  attack  from 
the  geometric  angle  of  incidence.  Nevertheless,  these  tests  are  often  used  in 
two-dimensional  comparisons  because  of  a  general  iack  of  expet  imental  data,  and 
for  that  reason,  we  do  so  also. 
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Comparison  of  pressure  distributions  on  the  NACA  4412  section  at  an  equivalent 
tv/o-dimensional  flow  angle  of  6.4  degrees  (geometric  incidence  of  8  degrees)  is 
shown  in  Figure  12,  Agreement  on  the  upper  surface  is  excellent,  but  the  lower 
surface  shows  slight  differences.  As  expected,  there  is  disagreement  near  the 
trailing  edge  caused  by  the  thickening  boundary  layer.  Also  shown  in  Figure  12  is 
the  pressure  distribution  computed  from  ordinates  interpolated  (in  the  pressure  program) 
from  the  tabulated  measured^  offsets  of  the  foil.  Some  of  the  humps  and  hollows  in  the 
measured  pressure  distribution  are  better  predicted,  but  a  new  predicted  hump  is  also 
obtained  at  quarter  chord.  (The  measured  nondimensional  ordinate  at  this  point  is 
0.0012  greater  than  the  computed  value.)  Essentially  though,  the  pressure  distribution 
is  the  same  as  for  the  mathematical  ordinates. 

Reference  3  contains  one  of  the  few  tabulated  sets  of  measured  two-dimensional 
pressure  distributions  with  a  sufficient  number  of  experimental  points  near  the  nose. 

The  model  tested  was  a  symmetric  RAE  101  section  with  a  30-inch  chord  and  10-percent 
thickness  ratio.  The  large  model  (surfaces  were  accurate  to  -t  0.0003  c)  and 
accurately  measured  angle  of  incidence  (to  the  nearest  0.01  degree)  makes  these 
tests  valuable  for  comparative  purposes. 

Two  Reynolds  numbers  were  considered  in  the  experiments.  For  comparison, 
we  have  selected  the  lower  (R  =  1 .6  x  ’0® )  since  it  corresponds  to  a  lower  Mach 
number  (U  =  100  ft/sec).  (It  is  not  completely  clear  whether  or  not  compressibility 
corrections  were  applied  to  the  data.)  Measured  and  predicted  pressure  distributions 
are  shown  in  Figures  13,  14,  and  15  for  three  angles  of  attack:  0,  2.05,  and  4.09 
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Figure  i3  —  Measured  and  Predicted  Pressures  on  the  RAE  101  Section,  10-Percent 

Thick.  Cl  -  0.  a  -  0 


Figure  14  —  Measured  and  Predicted  Pressures  on  the  RAE  101  Section,  10-Percent 

Thick,  Cl  =  0.218,  a=  2.05  Degrees 


Figure  15  ~  Measured  and  Predicted  Pressures  on  the  R.AE  101  Section,  10-Percent 

Thick,  Cl  =  0.430,  a  =  4.09  Degrees 
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degrees.  Agreement  is  excellent  except  at  the  trailing  edge.  This  deviation  is 

again  due  to  the  rapid  thickening  of  the  boundary  layer  at  the  trailing  edge. 

COMPARISON  WITH  MEASURED  CAVITATION  INCEPTION 

Although  measured  pressure  distributions  are  hard  to  find,  incipient  cavitation 

tests  have  been  conducted  in  several  laboratories.  With  the  classical  assumptions 

that  cavitation  starts  when  the  local  pressure  falls  to  the  vapor  pressure  of  the  flowing 

liquid,  cavitation  inception  may  be  predicted  from  the  pressure  distribution  since  the 

P  “  P 

cavitation  index  a  - -  is  equal  to  minus  the  minimum  pressure  coefficient 

'  \/2p\P  ^ 

P  ■  P  • 

X  rn  I  n 

■•Cp  =  — — — — —  .  It  is  important  to  note  that  unsteady  effects  are  not  considered 
when  computing  the  minimum  pressure.  That  'S,  the  testing  is  actually  done  for  either 


varying  P^^or'i^"’  whereas  the  predictions  are  for  fixed  flow  conditions.  This  means 


that  there  's  some  doubt  that  the  predicted  miminum  coincides  with  the  actual  test 
minimum.  Assuming  that  they  are  equal  is  the  usual  "quasi-steady"  approach,  generally 
the  only  reasonable  bar'Ts  of  solution.  For  the  (assumed)  slow  variations  in  flow 
conditions  here,  it  is  reasonable  to  expect  they  would  be  equal. 

Cavitation  inception  data  from  two  different  laboratories  will  be  considered: 

(1)  tests  by  the  California  Institute  of  Technology  (CIT)^’  on  the  NACA  4412  section 
for  various  Reynolds  numbers  and  (2)  tests  by  Vosper  Limited^”^  on  elliptic-parabolic 
sections  of  various  thickness  and  camber  ratios  (NACA  a  =  1  camberline)  at  constant 
Reynolds  number  {R^  =  1x10  ). 

Measured  values"^  of  q-.  for  (visual)  inception  and  disappearance  of  covitation 


at  two  different  Reynolds  numbers  are  compared  with  measured^  volues  of  -C 


P  . 

mm 
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for  the  NACA  4412  in  Figures  16a  and  16b.  Predicted  minimum  pressures  wore 

computed  for  values  of  V  and  o  interpolated  from  subcavitating  tests  of  the 

°e 

foil  used  in  the  'nception  tests. 

The  good  agreement  of  predicted  and  measured  -Cp  is  surprising  since  the 

^min 

angles  of  zero  lift  and  lift-curve  slopes  are  not  the  same  in  the  various  tests.  Also, 
for  negative  angles  of  attack,  the  measured  pressure  distribution  was  not  token  at  a 
sufficient  number  of  points  to  ensure  obtaining  the  minimum  pressure. 

Cavitation  prediction  from  computed  minimum  pressure  is  not  in  such  good  agreement 
(the  pfedictiun  is  conservative)  although  agreement  improves  with  increasing  test 
free-stream  velocity.  Some  of  the  discrepancy  may  be  attributed  to  nonstecdy  effects 
(resulting  in  a  poor  prediction  of  the  minimum  pressure)  although,  as  already 
discussed,  such  effects  should  be  small.  Again  the  fault  seems  to  be  with  the  assumption 
that  cavitation  begin;'  when  the  local  pressure  falls  to  the  vapor  pressure.'*'  Perhaps 
equally  inpcrtanr  in  accounting  for  the  inaccurate  prediction  is  the  small  size  of  the 
mode!  (3-irich  chord).  For  such  a  small  triodel,  slight  machining  errors  could  result 
in  iatge  changes  from  the  computed  pressure  distribution.  It  would  thus  be  of  con¬ 
siderable  interest  tc  have  measured  ordinates  rather  than  ordinates  for  only  the 
mathematical  foil. 

Machining  should  be  most  inaccurate  for  the  rapidly  changing  geometry  of  the 
nose,  and  at  negative  angles  of  attack,  the  minimum  pressure  is  close  to  the  nose. 

■*'This  hos  been  known  for  a  considerable  time  although  it  is  generally  ossumed  that 
for  "engineering  problems,"  conditions  are  such  that  the  assumption  Is  adequate. 

(See  References  23  and  24.) 
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Figure  16  -  Comparison  of  Incipient  Cavitation  Number  and  Minimum  Pressure  Coefficient, 

NACA  4412 
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It  is  possible,  then,  that  machining  inaccuracies  also  contribute  to  the  large  differences 
at  negative  incidence. 

Since  the  incipient  cavitation  number  for  two-dimensional  hydrofoils  increases 
with  increasing  free-stream  velocity  and  since,  in  “certain  instances, "  increasing 
the  body  size  also  increases  the  incipient  cavitation  number,^®  it  is  imperative  that, 
with  cur  present  lack  of  knowledge  of  the  proper  scaling  lav/s,  laboratory  testing 
should  be  done  for  environments  approaching  full-scale  conditions  if  meaningful  data 
are  to  be  taken,  A  better  alternative  is  to  develop  realistic  scaling  laws  so  that 
models  of  small  size  may  be  tested  at  low  Reynolds  numbers. 

The  cambered  9-percent  thick  foil  was  selected  for  comparison  with  the  Vosper 
test  data.  In  these  tests,  only  (visual)  cavitation  disappearance  was  recorded. 
Experimental  and  predicted  values  are  shown  in  Figure  17.  Figure  17a  gives  the 
raw  data  points  and  Figure  17b  the  faired  data  corrected  (by  Vosper)  for  tunnel 
interference.  The  raw  data  are  included  to  show  the  experimental  scatter.  In  both 
cases,  the  prediction  is  based  on  the  experimental  lift-curve  slope  corrected  for 
tunnel  interference.  The  raw  data  show  that  the  agreement  for  angles  greater  than 
zero  is  generally  good  except  for  the  one  point  at  1/2  degree.  For  negative  angles 
of  attack,  the  prediction  is  again  too  conservative. 

Although  the  model  chord  was  relatively  large  (8  inches)  in  these  tests,  the 
machining  was  accurate  to  only  ±  0.005  inch®^  or  ±  0.0006c.  From  actual  numerical 
test  cases  run  using  the  pressure  program,  fiiindom  differences  of  this  size  in  the 
ordinates  could  account  for  the  unusual  result  that  a  >-  C^,  near  zero  angle 
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at  small  incia'ences  was  also  found  for  other  foils 


of  attack. 


(However, 


rT;>- 


mm 

tested  by  Vesper.  It  is  interesting  to  noie  that  the  air  content  during  these  tests  was 
only  50  percent  of  saturations^  whereos  air  content  for  the  CIT  tests  was  about 
35  percent  of  saturation. s-) 

Although  measured  ordinates  would  be  useful  in  determining  more  precisely  the 

differences  betv/een  q.  and  -  Cn  ,  the  differences  are  too  great  to  be  attributed 

‘  min 

solely  to  inaccurate  machining.  Most  of  the  differences  seem  o  result  of  assuming 
that  cavitation  occurs  when  the  local  pressure  falls  to  the  vapor  pressure. 


SUMMARY  AND  CONCLUSIONS 


A  numerical  method  for  calculating  the  two-dimensionol  pressure  distribution 
on  arbitrary  profiles  with  arbitrary  lift  has  been  explained.  The  development  is  based 
on  an  empirically  modified,  approximate,  conformal  transformation  of  the  circle  and 
is  limited  to  flow  conditions  before  stall.  In  all  cases  considered  so  for,  the  numerical 
approximate  method  gives  integrated  lift  coefficients  that  differ  by  less  than  2  percent 
from  the  assigned  value. 

For  the  examples  tested,  when  the  pressure  distribution  about  a  foil  at  a  given 
incidence  is  computed  with  the  appropriate  lift,  agreement  is  good  with  both  exact 
potential  theory  and  experimental  results.  Comparisons  with  measured  pressure 
distributions  show  considerable  disagreement  near  the  trailing  edge  caused  by  the  thicken¬ 
ing  boundary  layer.  The  assumed  empirical  modification  introduced  in  the  potential 
theory  cannot  adequately  represent  the  flow  conditions  near  the  trailing  edge  since 
it  assumes  a  stagnation  point  there  in  contrast  with  the  experimental  result  of  almost 
free-stream  pressure. 


Although  predicted  pressure  distributions  agree  well  with  measurement,  comparisons 
of  the  minimum  pressure  coefficient  and  the  incipient  cavitation  index  show  differences. 
In  the  example  considered,  the  differences  decrease  with  increasing  test  Reynolds  number 
(i.e.,  rree-stream  velocity).  Also,  covitation  prediction  on  the  upper  surface  is  better 
than  on  the  lower  and  is  generally  conservative  on  both.  Some  of  the  discrepancy 
may  be  explained  by  machining  inaccuracies,  especially  near  the  nose.  Hence,  it  would 
be  of  interest  to  have  measured  ordinates  to  use  in  the  pressure  program  when  making 
comparisons  with  the  test  results.  If  we  assume  that  nonsteady  effects  in  the  testing 
were  small,  most  of  the  differences  seem  a  result  of  cavitation  not  occuring  when 
the  minimum  pressure  is  equal  to  the  vapor  pressure.  This  points  out  the  need  for  a 
better  scaling  law  if  laboratory  data  are  to  be  useful  in  predicting  full-scale  results. 

This  investigation  has  called  attention  to  two  areas  of  deficiency  in  experimental 
results.  First,  there  is  a  lack  of  pressure  distributions  on  two-dimensional  cambered 
foils.  Second,  cavitation  inception  tests  have  not  been  carried  to  the  point  of 
constant  (j.  for  increasing  Reynolds  number.  Further  investigation  in  both  of  these 
areas  would  be  valuoble  to  further  confirm  the  calculation  method  of  this  paper. 
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APPENDIX  A 


EXAMPLES  OF  EXACT  CONFORMAL  TRANSFORMATION 
If  specific  numbers  are  substituted  for  the  transformation  coefficients  in  Equations 
[2jand  [31,  we  will  obtain  foil  shopes  from  exact  theory  which  can  be  used  to  check 
the  approximate  method  where  x  has  been  approximated  by  Equation  [41 
SYMMETRICAL  FOILS 

Consider  first  symmetrical  shapes.  These  forms  will  be  generated  if  all  are 
zero.  For  our  purposes,  sufficient  generality  is  obtained  by  setting  ail  A^,  n  >  2, 
to  zero.  Then  we  have 

X  =  A^+(A_|+A^)  CoS  0  +  A  2  cos  2^ 

Y  “  (A  ^  -  A  ^  )  sin^  -  A  ^  sin  2tp 

Letting  Y  =  c  (sin  0  -  5  sin  20)  and  requiring  x  to  lie  between  0  and  1 ,  v/e  obtain 

X  =  1/2  (1  +  cos^)  +  e6  {cos  20-  1) 

From  Equation  04],  the  exact  velocity  is  obtained: 

[  sm^  cosa+  (1  -  cos0)  sinal 

Also,  from  Equation  [12] the  lift  coefficient  is 

C^  =  2tx  (  1  +  2  e)  sina 
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First,  consider  the  case  6=0,  then 


Y  =  esin 

X  =  1/2(1+  eoso) 

which  is  the  ellipse  of  thickness  ratio  r  =  2  e  . 

The  velocity  and  lift  coefficient  are 

„  sin</)cosa  +  (1  -  coso)  sina 

-a-=  (1+t)  ,  ,  z: - ^ 

U  y  sin^</)+'^  cos^^p 

=  2m  (  1  +  x)  sin  a 

(Note  that  this  includes  the  special  case  of  the  flat  plate,  r  =  0.) 

Since  x  is  precisely  the  form  of  the  approximation,  the  velocity  from  the 
approximate  method  is  also  exact. 

For  the  case  of  nonzero  e  and  6,  we  can  find  (P  in  terms  of  x: 

V'/4-3  6e0/2-2e6-x)  -  1/2 

(/)  =  arc  cos 

4  6  e 

and  obtain  ^  values  which  give  x,  Y,  and  q/U  at  the  x  values  used  in  the  numerical 
method.  This  has  been  done  for  four  cases: 


6  =  1/2 


0.2 

0.1 


6  =  1/4 


jT=  0.2 

I  T  =  0.1 


For  6  =  1/2,  the  foil  has  a  cusped  tail  similar  to  the  Joukowski  foil.  Maximum 


Y  value  occurs  at  ^  -  120  degrees,  from  which  e 


3t/T 


.  A  comparison  of  exact 


and  approximate  values  is  given  in  Part  (A)  of  the  following  table. 


56 


TABLE  A-1 

Comparison  of  Exact  and  Approximate  Velocity  Distribution  for  Two  Foils 


A)  For  the  Foil: 


X  =  4  (t  »  cos  o)  *  — (cos  20  -  1) 
2  3/3 

Y  =  (sm  0  -  i  sin  201 
3/3  2 


T  =  .20 

T^.IO  1 

X 

Y_ 

* 

c 

exact 

c 

''approx 

®exact 

n 

“^approx 

MM 

c 

exact 

c 

approx 

D 

exact 

D 

approx 

0 

0 

0 

■I 

;.4952 

6.4036 

0 

0 

0 

13.9904 

12.9480 

.007596 

.031655 

.7370 

mm 

7.0136 

5.9955 

.014381 

.9340 

.9343 

9.8212 

9.0912 

.030154 

.059631 

1.2296 

5.8179 

4.9847 

.027447 

1.1659 

1.1662 

6.0849 

5.6346 

.066987 

.081097 

1.4283 

mm 

4.4618 

3.8352 

.038080 

1.2116 

1.2119 

4.1635 

3.8577 

.116978 

.094556 

1.4422 

1.4454 

3.3303 

2.8741 

.045505 

1.2072 

1.2075 

3.0565 

2.8341 

.178606 

.099855 

1.3793 

1.3828 

2.4984 

2.1657 

.049373 

1.1846 

1.1849 

2.3431 

2. 1746 

.25 

.097881 

1.2979 

1.3004 

1.9067 

1.6603 

.049755 

1.1547 

1.1549 

1.8465 

1.7154 

.328990 

.090159 

1.2182 

J..2204 

1.4824 

1.2965 

.047083 

1.1225 

1.1228 

1.4816 

1.3779 

.413176 

.078485 

1.1479 

1.1499 

1.1708 

1.0282 

.042038 

1.0908 

1.0910 

1.2027 

1.1197 

.5 

.064652 

1.0882 

1.0900 

.9351 

.8244 

.035440 

1.0609 

1.0612 

.9826 

.9157 

.586824 

.050267 

1.0386 

1.0402 

.7516 

.6649 

.028131 

1.0339 

1.0341 

.8043 

.7502 

.671010 

.036642 

.9977 

.9992 

.6045 

.5364 

.020881 

1.0099 

1.0101 

.5562 

.6127 

.75 

.024738 

.9645 

.9659 

.4832 

.4299 

.014317 

.9893 

.9894 

.5304 

.4956 

.82»394 

.015149 

.9379 

.9392 

.3804 

.3392 

.008881 

.9719 

.9721 

.4212 

.3939 

.883022 

.008111 

.9171 

.9184 

.2909 

.2598 

.004805 

.9579 

.9581 

.3242 

.3033 

.933013 

.033540 

.9016 

.9028 

.2108 

.1886 

.002114 

.9471 

.9472 

.2361 

.2210 

.969846 

.001074 

.0908 

.8919 

.1372 

.1228 

.000645 

.9394 

.9396 

.1542 

.1443 

.992404 

.000136 

.8844 

.8855 

.0676 

.0506 

.000082 

.9348 

.0713 

1.0 

0 

0 

0 

0 

0 

0 

°  . 

0 

B) 


For  the  Foil: 


X  =  4  (1  cos  o)  4"  (cos  2o  -  1) 

i  4 

Y  -  €(sm  0  -  -  sin  2o) 

4 


e  -  .45417  r 


T  .20 

T  =  .10 

X 

Yt 

c 

exact 

c 

approx 

^exact 

^approx 

Yt 

("exact 

("approx 

(*exact 

(^approx 

0 

0 

0 

0 

8.6728 

7.9168 

0 

0 

0 

16.0122 

15.3005 

.007596 

.025944 

.7397 

.7403 

7.6502 

6.9847 

.012328 

.9377 

.9378 

10.2197 

9.7659 

.030154 

.049890 

1.1375 

1.1384 

5.8403 

5.3350 

.023823 

1.1206 

1.1207 

6.0605 

5.7921 

.066987 

.070132 

1.2895 

1.2904 

4.3603 

3.9865 

.033746 

1.1585 

1.1586 

4.1239 

3.9420 

.116978 

.085472 

1.3270 

1.3280 

3.3070 

3.0269 

.041532 

1.1621 

1.1622 

3.0462 

2.9125 

. 178606 

.095321 

1.3143 

1.3152 

2.5600 

2.3461 

.046838 

1.1532 

1.1534 

2.3602 

2.2572 

.25 

.099687 

1.2803 

1.2812 

2.0171 

1.8513 

.049562 

1.1387 

1.1388 

1.8827 

1.8012 

.328990 

.099068 

1.2388 

1.2397 

1.6117 

1.4813 

.049827 

1.1214 

1.1215 

1.5294 

1.4637 

.413176 

.094293 

1.1964 

1.1972 

1.3008 

1.1973 

.047942 

1.1032 

1.1033 

1.2560 

1.2024 

.5 

.086389 

1.1563 

1.1571 

1.0565 

.9738 

.044342 

1.0850 

1.0851 

1.0369 

.9931 

.586824 

.076410 

1.1201 

1.1208 

.8598 

.7936 

.039526 

1.0678 

1.0678 

.8565 

.8205 

.671010 

.065341 

1.0881 

1.0888 

.6979 

.6449 

.033997 

1.0518 

1.0519 

.7043 

.6750 

.75 

.054009 

1.0605 

1.0611 

.5615 

.5194 

.028205 

1.0375 

1.0375 

.5730 

.5493 

.821394 

.043034 

1.0369 

1.0375 

.4438 

.4109 

.022512 

1.0249 

1.0250 

.4573 

.4385 

.883022 

.032809 

1.0166 

1.0171 

.3399 

.3149 

.017162 

1.0141 

1.0142 

.3533 

.3388 

.933013 

.9976 

.9982 

.2457 

.2278 

.012279 

1,0047 

1.0048 

.2577 

.2472 

.969846 

.9733 

.9739 

.1578 

.1464 

.007867 

.9'J48 

.9949 

.1679 

.16U 

.992404 

.007356 

.8984 

.8988 

.0723 

.0670 

.003830 

.9696 

.0812 

.0779 

1.0 

0 

0 

0 

0 

0 

0 

0 

0 

0 
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For  6  =  1/4,  the  trailing  edge  is  rounded  and  the  general  foil  shape  is  similar 

is-]  ^ 

to  conventional  foils.  The  maximum  value  of  Y  occurs  at  (p=  (tt  -  arc  cos - ), 

2 

for  which  e  =  0.45417  r .  A  comparison  with  the  approximate  solution  is  given  in 
Part  (B)  of  the  tab’e. 

These  examples  show  negligible  differences  in  the  C  {(p)  component  of  velocity 
for  the  10-percent  thick  foils.  For  the  20-percent  thickness  forms,  the  differences 
are  larger  but  still  slight.  It  is  interesting  to  note  that  the  approximate  velocity  is 
greater  than  the  exact  in  this  case.  The  D((D)  component  of  velocity  shows  large 
differences,  but  the  approximate  velocity  is  lower  than  the  exact  in  this  case,  which 
is  the  trend  indicated  by  viscous  effects.  Also  this  is  the  component  of  velocity 
which  is  multiplied  by  sin  a  so  that  its  contribution  to  the  total  velocity  is  small 
in  the  normal  incidence  range. 

In  general,  these  comparisons  indicated  that  there  are  larger  errors  in  the 
approximate  computations  for  foils  which  show  large  departures  from  an  ellipse. 

Most  foils  in  use  today  have  their  maximum  thickness  near  mjdchord  and  do  not  have 
cusped  tails;  both  of  these  properties  contribute  to  accurate  calculations. 

CAMBERED  FOILS 

For  foils  with  camber,  we  consider  terms  in  Equation  [21  and  fSl  to  the 
second  harmonic,  i.e., 

X  =  A  +  (A  +  A  )  cos  o  +  (B  -  B  )  sirjtp  +  A  cos  2fp  +  B  sin  ^ 
0-11  1-1  2  '  2 

Y  =  B  +  (A  -  A  )  sin  y)  +  (B  +  B  )  cosP  -  A  sin  2^  +  B_  cos  2(p 
O  - 1  I  - 1  I  2  2 
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We  take 


Y  =  e  (sinO-  5  sin  2(pj  -f  —  (  1  -  cos  2o)  +  y  (1  -  cosip) 

2 

dx  \ 

and  set - j  =  0,  x(0)  =  1 ,  x(tt)  =  0  to  obtain 

dcp  !  fp  =  Q 

x  =  1/2(1+  cosp)  +  e  6 (cos  2tp-  1)  +  f  sin^  -  —  sin  2y> 

dx  ' 

Unfortunate! y,’'  these  foils  do  not  have -  =  0  at  the  leading  edge  which  means 

dfp 

there  is  some  overhang  there.  The  constant 7  is  included  to  raise  the  nose  to  lie 
close  to  the  x  axis.  The  velocity  distribution  on  these  foils  is 


iL  - 
U  " 


cos  a  c)  sinf/5+  (f +  y)  (1  *■  COS0)  ] 

+  sing  [(1/2  +  c)  (1  -  cos^)  -  (f +  y)  sinrp  ] 


dx  dy  ® 
dp  '  dy) 


dx 


where - "•  -  1/2  sin^  +  f  cos'p  -  2  e5  sin  2y)  -  f  cos  2^ 

ckp 


dy 

— ^  ?  (cosy)  -  26  cos  2<p)  +y  siny)+  f  sin  2p 

dup 


1  +  2  c 

The  lift  coefficient  Is  C  =  2-n  — 'z 

L  ^ 


1  + 


X  =  -  arctan 


2(f  +  y) 

1  +  2  e 


2(f  +  y) 
1  +  2  e 


sin  (a-  x)/ 


The  complexity  of  the  x  {(p)  term  does  not  permit  easy  inversion  to  obtain 
‘P  =  y5(x).  Instead,  we  simply  calculated  x,  Y,  and  q/U  at  many  rp  values  and  put 


*Nor  do  these  expressions  permit  us  to  put  e  =  0  ond  obtain  simply  a  camberline 
shape  since  the  x  relation  would  then  give  a  curve  which  crosses  itself  and  hence 
is  meaningless. 
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fhe  (x,Y)  relation  into  the  arbitrary  input  section  of  the  pressure  prc^ram.  Since  we 
were  unable  to  calculate  the  exact  velocity  at  the  points  of  fhe  numerical  method, 
the  comparisons  must  be  made  graphically. 

2 

The  thickness  portion  of  the  foil  was  taken  as  Yp  = — (sine?  -  1/2  sin  2<o) 

3y3 

with  T  =  0. 10.  The  constants  f  and  7  were  first  taken  as  0.03  and  0.01/  respectively, 
which  gives  a  camber  ratio  of  about  0.0369  and  then  f  and  7  were  taken  as  0.015  and 
0,0075  which  gives  a  camber  ratio  of  obout  0.0185.  The  foil  shapes — rotated  end 
shrunk  to  put  the  nose  at  (0,0) — and  velocity  distribution  computed  from  exact  theory 
and  the  approximate  theory  are  shown  in  the  following  figures  for  both  cases  at  two 
angles:  0  and  5  degrees  (angles  referenced  to  the  unrotated  foil  shape).  The 
calculations  show  the  approximate  m.ethod  is  sufficiently  accurate  for  practical  use 
at  0  degrees.  At  5  degrees,  the  opproximation  is  showing  enough  error  to  be 
questionable  from  a  potential  standpoint.  However,  the  error  is  ago  In  in  the  direction 
indicated  for  viscous  effects. 

Also  shown  In  these  figures  is  the  velocity  distribution  computed  with  A  a 
adjusted  to  give  the  same  lift  os  the  exact  theory.  These  comparisons  indicate  that  the 
approximate  method  is  sufficiently  accurate  for  most  work  when  the  lift  coefficient 
is  known. 
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NON-DIMENSIONAL  ORDINATE 


APPENDIX  B 


DISCUSSION  OF  AIRFOIL  GEOMETRY 

In  many  cases,  a  thickness  distribution  is  combined  with  a  camber!  ine  by  laying 
off  the  thickness  ordinate  perpendicular  to  the  camberline. One  result  of  this 
combination  method  is  a  nonzero  ordinate  at  the  leading  edge  (the  center  of  the  nose 
radius  lies  along  the  camberline  tangent  at  the  leading  edge^^.  A  taim  may  be 
added  to  Equation  {241  to  include  the  nonzero  ordinate  at  the  nose.  However,  when 
that  //as  done  for  the  NACA  4412,  the  computed  pressure  distributions  gave  a  greater 
negative  peak  for  positive  incidences  than  either  the  experimental  or  exact  potential 
values.  This  result  could  be  anticipated  from  c  consideration  of  the  approximations 
in  the  theory.  At  the  nose,  the  velocity  due  to  thickness  is  zero  and  the  only 
velocity  is  that  of  the  camberline.  Since  the  velocity  at  a  point  is  considerably 
influenced  by  the  ordinate  at  that  point,  a  large  ordinate  at  the  nose  means  a  large 
camber! ine  velocity  and  correspondingly  large  errors  in  the  computed  velocity.  To 
handle  this  nose  ordinate,  the  profile  can  be  rotated  to  put  the  nose  ct  (0,0).  When 
the  pressure  distribution  was  computed  for  c  rotated  foil,  agoin  the  4412  section, 
reasonable  results  were  obtained  although  the  predictions  for  potential  flow 
(  Aa  =  0)  were  somewhat  lower  than  the  exact  results  (see  Figure  11),  as  expected 
from  the  comparisons  in  Appendix  A. 

in  order  to  rotate  the  foil,  it  is  first  necessary  to  determine  the  nose  point 
accurately.  Points  along  the  profile  are  determined  from  the  expressions* 

*These  expressions  would  follow  more  logically  if  the  abscissa  for  the  thickness  ordinate 
were  measured  along  the  camberline.  Presently,  this  is  not  done.^® 
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upper  surface:  (x  -  Yy  sirs  S,  Y^t  Yj  cos  0) 

lower  surface:  (x  +  Y  sin  9,  Y  -  Y  cos  0)  [B-l] 

1  I 

where 


is  the  thickness  ordinate. 


is  the  camberline  ordinate,  and 


9  is  the  camberljne  inclination,  8=arctan 

However,  these  expressions  do  not  give  a  value  for  the  nose  point.  The  nose  can 
be  found  by  recalling  that  the  center  for  the  nose  radius  lies  along  the  camberline 
tangent  at  the  leading  edge,*  as  below 


Hence,  the  nose  (the  point  of  minimum  x)  is  given  by 


x  =  0 


X  —  -  P  (1  ~  cos  0  ) 

N  LE^ 


^LE^'"  ®LE 


[6-21 


The  nose  can  be  put  at  (0,0)  by  the  coordinate  system  translation  (maintaining  a 
chord  length  of  unity): 


*See  footnote  on  preceding  page. 
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-X 


X  =■ 


N 


]  -X 


N 


Y-  Y 


Y  = 


N 


1  -X 


N 


Then  the  coordinates  con  be  rotated  through  the  angle 

Y. 


U)  =  arctan 


N 


1  -X 


N 


[B-3'» 


to  measure  the  ordinates  from  the  nose-tail  line.  The  new  coordinates  with  a 
chordlength  of  unity  are 


X  cos  r£  -  Y  sin 


[B-5] 


In  the  pressure  program,  the  rotated  ordinates  are  found  from  Equation  IB-Sl. 

Although  the  above  expressions  are  exact,  they  do  require  interpolation  between 
the  computed  values,  or  an  iteration,  to  find  the  ordinates  for  a  fixed  station. 

An  alternative  is  an  approximate  equation  for  combining  a  camberline  and  thickness 
distribution  for  fixed  x.  Such  an  expression  may  be  obtained  by  expanding  the 
difference  between  the  known  ordinate  in  Equation  IB-1]  and  the  unknown  ordinate 
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at  X  in  a  Taylor  series.  The  results,  retaining  only  the  lowest  order  terms,  ore 


V!  +Y^'  Y^' 


Y  =  Y  -  Y 
L  C  T 


-Yt'Yc' 


rB-61 


where  primes  denote  differentiation  with  respect  to  x. 

With  Y^  =  2  Vp  X  +  . . .  these  expressions  give  the  nose  as  {0,  Yj^)  where 

Yn='’leYc'«»- 


The  rctorion  to  put  the  nose  at  (0,0)  is  approximately 


Y|j«Y-Y^(l -x)  rB-71 

These  expressions  are  quite  accurate  for  small  thickness  and  camber  ratios. 

However,  they  are  not  used  in  the  pressure  program  and  are  included  only  to  show 
a  method  of  combining  a  thickness  distribution  perpendicular  to  a  camberline  at 
fixed  X  so  that  ordinates  may  be  obtained  at  the  stations  required  in  the  pressure 
program.  As  already  mentioned,  there  is  a  considerable  saving  of  machine  time  when 
the  ordinates  are  given  at  the  required  stations. 


65 


APPENDIX  C 


LISTING  OF  THE  FORTRAN  STATEMENTS 
FOR  THE  PRESSURE  PROGRAM 
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